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Abstract 

The system of compressible adiabatic flow through porous media is consid- 
ered in M 3 in the present paper. The global existence and uniqueness of classical 
solutions are obtained when the initial data is near its equilibrium. We also show 
that the pressure of the system converges to its equilibrium state at the same 
L 2 -rate (1 + i)~ 3 as the Navier-Stokes equations without heat conductivity, but 
the velocity of the system decays at the L 2 -rate (1 + i)~ 3 , which is faster than 
the L 2 -rate (1 + for the Navier-Stokes equations without heat conductivity 
[3]. 

Mathematics Subject Classification (2000). 76W05, 35Q35, 35D05, 76X05. 

Keywords. Euler equation with damping, Global existence, Large time behav- 
ior. 



1 Introduction 

The motion of compressible adiabatic flow can be modeled by the compressible Euler 
equations with frictional damping terms, say, the following balance laws: 

d t p + V -(pu) = 0, 

d t (pu) + V • (pu <g> u) + Vp = —apu, (1.1) 
dt{pS) + V • {pu£ + up) = -apu 2 , 

for (t,x) G [0, +<x>) x M 3 . Here p, u = (iti, 1*2, ^3)* and p denote the density, the 
velocity, and the pressure respectively. The total energy £ = ^- + e, where e is the 
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internal energy. The constant a > models friction. In this paper, we will consider 
only polytropic fluids, so that the equations of state for the fluid is given by 

p = Rp9, e = c v 9, 

where R > 0, c v > are the universal gas constant and the specie heat at constant 
volume respectively. 

In the case of isentropic flow where s=const., (1.1) takes the form 

f d t p + V • (pu) = 0, ^ _ 

[ d t (pu) + V • {pu <g) u) + Vp = —apu. 

For one dimension, the system (1.2) is hyperbolic with two characteristic speeds A = 
u ± a/p'(p)- As a vacuum appears, it fails to be strict hyperbolic. Thus, the system 
involves three mechanisms: nonlinear convection, lower-order dissipation of damping 
and the resonance due to vacuum. The global existence of a smooth solution with small 
data was first proved by Nishida [23,24], and the behavior of the smooth solution was 
studied in many papers, the reader can refer to [1,2,5-7,10,12-15,18,20,21,30,32] and 
references therein. In mult i- dimension space, Wang- Yang [29] considered the time- 
asymptotic behavior of solutions in multi-dimensions, the global existence and point- 
wise estimates of the solutions were obtained. Sideris-Thomases-Wang [27] proved that 
damping prevents the development of singularities in small amplitude classical solutions 
in three-dimensional space, using an equivalent reformulation of the Cauchy problem 
to obtain effective energy estimates. Fang-Xu [4] studied the existence and asymptotic 
behavior of C l solutions to the multi-dimensional space on the framework of Besov 
space. Jang-Masmoudi [16] studied well-posedness of compressible Euler equations in 
a physical vacuum. The optimal estimates was obtained by Tan-Wu [28]. For initial 
boundary value problem, refer for instance to [10,26] and references therein. 

For the adiabatic flow where s ^ const., much less is known. For one dimension, 
the global existence of smooth solutions to the Cauchy problem has been proved in [11] 
and [31] for small initial data. The problem of large time behavior of these solutions 
is known only for some particular initial data; see [8,22]. For initial boundary value 
problem, refer for instance to [9,25] and references therein. 

From the physical point of view, the mult i- dimensional case of model (1.1) describes 
more realistic phenomena. Also the multi-dimensional compressible Euler equations 
carry some unique features, such as the effect of vorticity, which are totally absent in 
the one dimensional case and make the problem more mathematical challenge, system 
(1.1) and its time-asymptotic behavior are of great importance and are much less 
understood than its one dimensional companion. To our knowledge, there is not any 
work for the full system (1.1) in M 3 . 

The study of this paper is motivated by Duan-Ma [3], where the authors considered 
the compressible Navier-Stokes equations without heat conductivity in terms of the 
variables p, u and s. It is well-known that all thermodynamics variables p,9,e,p as 
well as the entropy s can be denoted by functions of any two of them. We take the two 
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variables to be p and s. Then the equation of state for the gas is then given by 

p = kp^+Re^^^, (1.3) 

where k > is a constant. Under the aforementioned assumptions, the system (1.1) in 
terms of the variables p, u and s reads 

d t p + -^J^pV • u + u ■ Vp = 0, 

d t u + (u ■ V)u + ^ = -au, (1.4) 
d t s + (u • V)s = 0, 

where p = p{p,s) is defined by (1.3). Notice that (1.4) is a hyperbolic system, where 
the dissipation comes from damping. We consider the initial value problem to (1.3) in 
the whole space R 3 with the initial data 

(p(x),u(x),s(x))\ t=0 = (po(x),U (x),S (x)) ► (poo,0,Soo) as \x\ — > oo, (1.5) 

where p^ > and Soo are given constants. 

Before we state the main results, let us introduce some notations for the use through- 
out this paper. C denotes some positive constant. The norms in the Sobolev Spaces 
if m (R 3 ) and W m ' q (M. 3 ) are denoted respectively by || • || m and || • \ \ miq for m > and 
q > 1. In particular, for m = we will simply use || • || and || • \ \lq- Moreover, we use 
(■, ■) to denote the inner product in L 2 (R 3 ). Finally, 

V = (d 1 ,d 2 ,d 3 ), di = d Xi , i = 1,2,3, 

and for any integer / > 0, V 1 f denotes all derivatives of order / of the function /. And 
for multi-indices a and j3 

a = (ai, a 2 , a 3 ), /3 = (/?i, /3 2 , /3 3 ), 

we use 

3 

i=i 

and = p^lpy where (3 < a. 

For the global existence and large time behavior of classical solutions, we have the 
following: 

Theorem 1.1. Let the initial data (po,v>o,so) be such that ||(po -Poo,«o, s o — s oo)||3 
is sufficiently small and \ \(po — Poo,U'o)\\l 1 is bounded. Then the initial value problem 
(1.4)-(1.5) admits a unique solution (p, u, s) globally in time with p > 0, satisfying 

p - Poo , u, s - Soo G C°(0, oo; H 3 (R 3 )) n C\0, oo; H 2 (R 3 )). 
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Moreover, there exists a constant C such that for any t > 0, 

-Poo,«)(t)||!+ Al|Vp(r)||| + ||«(r)||I)dr<C'o||(pd-Poo,tio)||I > 



||(s-Soo)(t)|| 3 <C'o||(po-Poo,«o,s-Soo)||3e C(,|l(po - p ° c ' ?t<,)11 ^^ 3 . 

Finally, there is a constant C\ such that for any t > 0, the solution (p, u, s) has the 
decay properties 

||(p-Poo)(*)||<Ci(l 

||V(p-p 0O )(t)|| 2 + ||w|] 3 < C!(l + *)-f, 
||a t (p,«,s)(t)||<C 1 (l + t)-i. 

The rest of the paper is organized as follows. We will reformulate the problem in 
Section 2. In section 3, we take Hodge decomposition to analyze linearized system and 
establish the L 2 time-decay rate. The proof of Theorem 1.1 is given in Section 4 and 
5. 



2 Reformulated system 

In this section, we first reformulate the problem as follows. Set 



K l = \ 7TT, s ) K 2 



(R + c v )p c 



where = p(poo, Soo)- Taking change of variables by 

(p, It, S) > (p + Poo, Kilt, S + Sqo), 

the initial value problem (1.4)-(1.5) is reformulated as 

d t p + k 2 V -u = F, 
d t u + k 2 Vp + aw = (7, 
<9 t s + ■ V)s = 0, 

(p,w, s)| t=0 := (po,Mo,So) ->■ (0,0,0) as |x| ->■ oo, 



(2.1) 



where 

F(p,u,s) = - (gigl^- pV -n-KiU-Vp, 
G(p,u,s) =-«i(S.V)«-^(i-^)Vp,. 

Here and in the sequel, for the notational simplicity, we still denote the reformulated 
variables by (p, u, s). 

Let us define the solution space and the solution norm of the initial value problem 
(2.1) by 

X(0,T) = {(p,u,s);p,u,se C°(0,T;# 3 (R 3 )) H C 1 ^, T; # 2 (M 3 )), 
Vp G L 2 (0, T; i^ 2 (IR 3 )), it G L 2 (0, T; # 3 (M 3 ))}, 



and 



N(0,T) 2 = sup \\(p,u,s)\\l+ [ T (\\Vp(t)\\l+\\u(t)\\l)dt, 

0<t<T Jo 
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for any < T < oo. The local existence and uniqueness of H solution can be estab- 
lished by following the methods in Kato [17] or Majda [19]. 

Proposition 2.1 (local existence). Let (po,u 0: s ) G H 3 (R 3 ) be such that 

inf {p {x) +Poo} > 0. 

Then there exists a positive constant To depending on N(0, 0) such that the initial 
value problem (2.1) has a unique solution (p, u, s) G X(0, T ) which satisfies N(0, T ) < 
2iV(0,0) and 

inf {p(t,x) +Poo} > 0. 

xGR 3 ,0<t<T 

To prove global existence of a smooth solution with small initial data, we establish 
global a priori estimates of the solution. 

Proposition 2.2 (A priori estimate). Let (po,u ,s ) e H 3 (R. 3 ) and (p ,u ) e L^R 3 ). 
Suppose that the initial value problem (2.1) has a solution (p,u,s) G X(0,T), where 
T is a positive constant. Then there exist a small constant e > and a constant C 2 , 
which are independent of T, such that if 

sup ||(p,u,s)(*)|| 3 <e, (2.2) 

0<t<T 

then for any t G [0, T], it holds that 

\\(p,u)(t)\\l+ A||Vp(r)||^ + |KT)||^r<C 2 ||(p ^o)||^ (2.3) 
Jo 

||^)||3<C 2 ||(po^o^o)||3e C2 ^, (2.4) 

where 

= ||(po,Uo)IUinff3- (2.5) 

Furthermore, there is a constant C3 such that for any t G [0,T], the solution (p,u,s) 
has the decay properties 

lb(*)ll < C7i(H-t)-4, (2.6) 
IIVKOII2 + ||«|| 3 < C^l + ^-t, (2.7) 
||a t (p,«,s)(t)||<Ci(l + t)-t. (2.8) 

Theorem 1.1 follows from Proposition 2.1 and Proposition 2.2 by standard conti- 
nuity argument. 
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3 Spectral analysis and linear L 2 estimates 



The linearized equations corresponding to system (2.1) is 
d t p + k 2 V -u — 0, 

d t u + K 2 Vp + au = 0, , , 

d t s = 0, {6A) 
(p,u,s)\ t=0 := (p ,u ,s ) ->■ (0,0,0) as \x\ ->■ oo. 

We take Hodge decomposition to analyze the first two of linearized system (3.1). 
For r G R, we denote by A r the pseudo differential operator defined by A r f = 
'^ r ~ 1 {\£,\ r f {£))■ Let v = A -1 div u be the "compressible part" of the velocity and 
oj = A~ 1 curl u (with [curl z)\ = djZ 1 — d.^) be the "incompressible part", then the 
first two of linearized system (3.1) writes 

d t p + k 2 Av = 0, 

d t v - n 2 Ap + av = 0, (3.2) 
d t w + aw = 0. 

Indeed, as the definition of v and u, and relation 

u = -A^Vv - A^div u (3.3) 

involve pseudo-differential operators of degree zero, the estimates in space H l (M. 3 ) for 
the original function u will be the same as for (v,ou). 

This section is devoted to the proof of the following results. 



Proposition 3.1. Let U = (po,u ) G H l (R 3 ) n L^IR 3 ), and U = (p,u) satisfies 
the first two equation of system (3.1). Let v := A~ l div u and uj := A _1 curl u. Then 
there exists a constant C such that for < \k\ < I, 

\\d k Mt)\\ < Ce- at \\U \\ k , (3.4) 

||^p(*)|| < C{\ + O-^^dl^olUi + Wok), (3.5) 
\\d k x v(t)\\ + \\d k x u(t)\\ 

< ca+^-t-^dic/oii^ + iic/oiu). 1 ' ' 

Proof: The estimate for oo is obvious, so let us focus on the first two equations of 
(3.2). In terms of the semigroup theory for evolutionary equation, the solution (p,v) 
of the first two equation of system (3.2) can be expressed via the Cauchy problem for 
V = (p, vf as 

V t = BV, V(0) = V o , t>0. 

which gives rise to 

V(t) = S(t)V =: e tB V , t > 0. 
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What is left is to analyze the differential operator B in terms of its Fourier expres- 
sion A and show the long time properties of the semigroup S(t). Taking the Fourier 
transform with respect to the space variable yields 

±V = m V with A(0 := ( -f! 

The characteristic polynomial of is A 2 + aX + k 2 |£| 2 and possesses two distinct 
roots: 

A ± (« = -|(l ± ^TMp!). (37) 
The semigroup e* A is expressed as 

where the project operators P± can be computed as 

P + = A( ?-y and P.= A «'-^. 

A_|_ — A_ A_ — A+ 

By a direct computation, we can verify the exact expression the Fourier transform 
G(£, t) of Green's function G(x, t) = e tB as 



G(£, t) =: e tA = e A +'P+ + e A -'P_ 

\+e x - t -\-e x + t K 2 |g|(e A + f -e A - t ) 



A+-A- A+-A- 
K 2 \(,\{e x + t -e x - t ) A + e A - t -A_e A + i _ q(e A + t -e A - t ) 
A-)- — A— A-)- — A— A-)- — A— 



(3.8) 



To derive the long-time decay rate of solutions whatever in I? framework or on 
point- wise estimate, we need to verify the approximation of the eigenvalues (3.7), we 
are able to obtain that it holds for |£| -C 1 that 



A+e A -* - A_e A +* -ae"^* + ^t e ~ at 



A+ — A_ —a 



e \ + t _ e \_t e -at _ e ^-t 



A i — A_ —a 



, (3.9) 



, (3.io) 



A + e A - f - A_e A+t a(e A +* - e A ~ f ) kIKI*-* „ 1 

A + -A_ A + -A_ ~ (1 "^ )e ' lel<<1 - (3 - U) 

To enclose the esimates, we also need to deal with the high frequency |£| ^> 1. By 
a direct computation, we have the following lemma: 
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Lemma 3.2. For any given constant rj > 0, there exist two positive constant C 
and i? which are just dependent of rj such that 

\\m,t)\\\<Ce-^ for |£| > 77, 

where 1 1 1 • 1 1 1 denotes the norm of matrix for simplicity. 

Let us now tackle the proof of (3.5) and (3.6). The proof relies on the use of explicit 
expression for p, v. By (3.8), we have 

A A+e A -* - A_e A +'„ /t 2 |£|(e A +' - e A -*) A 

p = a + - a_ Po w: — Uo - 

Combining (3.9), (3.10) and Lemma 3.2, we have the L 2 -decay rate for p as 

iipWHiw = Ji € |<„|p(e,*)i 2 de+ Ji € |>,ip(e,oi 2 de 

^2 it|2 



l|i«i 2 



(3.12) 



+Ce-^4 | ^(|p | 2 + |«o| 2 2 )^ 

+Ce-^\\{p ,u )\\l 2m 

< C(l + t)-i\\(p ,u )\\l 2mf]Lim . 

Taking the same argument, we have the L 2 -decay rate for v as 

Mt)\\h m = I K \< v \H^t)j 2 ^ + J^\^,t)\ 2 d^ 

+Ce-*°*f^(\p \* + \v \ 2 )dZ 

< C(l + t)-3||(p ,M )||i 2(R 3 )n£ i (R 3 ) . 

The L 2 -decay rate on the derivatives of (n, u) as 

\m>{t)\\h^) = i m ^mM^ +4>„ \z\ 2k \mt)?dt 

< Cj m<n U k e-^\\p,\ 2 + \h\ 2 )d^ (314) 

+Ce-^J^ v \e k (\Po\ 2 + \v m 

< C(l + t)-i- k \\(p ,u )\\ 2 Hkmr]Lim , 

and 

w&mhw = 4i<jei 2fc i^,t)i 2 ^+ii,i>,ier fc i^,oi 2 ^ 



^ Gf^l^e-^WW + ^m (3 15) 

+^ *^|>,le| 2fe (|Po| 2 + |^o| 2 )^ 
< C(l +t)-t- k \\(po,u )\\l 2mnLim , 
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for k > 1. By the relation between u and (v,oo), we can easily get the estimate for u: 

\\d k x u(t)\\ <C(l + t)-^ k \\U \\ HHR3)f]Ll 
The proof of proposition 3.1 is completed. 

We also need the following Sobolev's inequalities. 

Lemma 3.3. Let / G H 2 (R 3 ). Then it holds: 

(0 ll/ll^<C||V/||^||V/||F<C||V/|| i; 



(it) ||/|U 6 <C||V/||; 

(Hi) ll/IU^qi/Hi, 2<q<6. 

4 A priori estimates 

We suppose that the inequality (2.2) holds throughout this section and next section. 
The initial value problem (2.1) has a solution (p, u, s) in the space X(0,T) with some 
T G (0,+oo]. We also omit the variable t of all functions in the proof of different 
lemmas in this section for simplicity. 

In what follows, a series of lemmas on the energy estimates are given. Firstly, the 
energy estimate of lower order for (p, u) is obtained in the following lemma. 

Lemma 4.1. There exists a constant D x > suitably large which is independent 
of e such that 

^(^ik^,^)^)!! 2 + <^,^>(t)) + c(||v^(t)|| 2 + ||^(t)|| 2 ) < c||v^(t)|| 2 , (4.i) 

for any < t < T. 

Proof: Multiplying (2.1)!-(2.1) 2 by p, u respectively and then integrating them over 
M 3 , we have 

\ j t 1 1 (P, «) (t) 1 1 2 + a\ \u(t) 1 1 2 = (p, F) + (u, G) . (4.2) 

The two terms on the right hand side of the above equation can be estimated as follows. 
Firstly, for the first term, it holds that 



(p,F) = ~ {R+ ^ )K1 (p, pV -u) - «i (p, u ■ Vp) 
= (-i^hi + f )(p\v-u). 

It follows from Lemma 3.3, Holder's inequality and (2.2) that 

\(p 2 y-u)\ < ||p|| L 3||p|| L6 ||Vu|| kcWpUVpWWVuI 

< Ce(||Vp|| 2 + ||V«|| 2 ). 
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(4.3) 



(4.4) 



Putting (4.4) into (4.3,) we arrive at 

\(p,F)\ < C e (||Vp|| 2 + ||V M || 2 + |M| 2 ). (4.5) 
For the second term, we have 

\(u, G)\ < C(\(u, (u • V)«>| + \(u, (- - — )Vp>|). (4.6) 

P Poo 

Similar to the proof of (4.5), it follows from Lemma 3.3, Holder inequality and (2.2) 
that 

\(u,(u- V)u)| < ||«|| £ 3||u|| £ 6||Vu|| < C||u||i||Vu|| 2 < Ce\\Vu\\ 2 , (4.7) 

K«.(i-i)vp>l <IMlL.||J-illL.||Vp|| 

<C7||V«||||(p, a )||i||Vp|| (4.8) 
<Ce(||Vp|| 2 + ||V M || 2 ), 

where by (1.3) and (2.2), we have used the fact 

p~ Poo + 0(l)(p + s) 

and 

0(l)(p + s). 

P Poo 

Substituting (4.7) and (4.8) into (4.6), we obtain that the second term is bounded by 

\(u,G)\<Ce(\\Vp\\ 2 + \\Vu\\ 2 ). (4.9) 
Hence combining (4.2), (4.5) and (4.9) yields 

j t \\(p,u)\\ 2 + C\\u\\ 2 < Ce(||Vp|| 2 + ||V W || 2 ), (4.10) 

since e > is sufficiently small. 

Next we shall estimate ||Vp|| 2 . From (2.1) 2 , we have 

« 2 ||Vp|| 2 = (-ut, Vp) -a(u, Vp) + (G, Vp) . (4.11) 

By (2.1)i, the first term on the right hand side can be written as 



(-ut, Vp) = -f (Vp, u) + (Vp t , u) 

= -|(Vp,u)-(p t ,V-«) (4.12) 
= -i(Vp, if) + (k 2 V • u - F, V • u). 



It follows from the definition of F that 

|(-F,V-u)| < C(\(pV -u,V -u)\ + \(u- Vp,V-«)|) 

< C(||Vp||i||V-m|| 2 + ||Vp|| 1 ||Vm|| 2 + ||m||i||Vm|| 2 ) (4.13) 

< Ce||Vw|| 2 . 

10 



Taking the same argument to estimate (G, Vp), we have 

|(G,Vp)|<C e (||Vp|| 2 + ||V M || 2 ). (4.14) 
Using Cauchy-Schwarz inequality, we easily get 

\a(u,V P )\ < C\\u\\ 2 + 'jWVpW 2 . (4.15) 

Since e > is small enough, putting (4.12), (4.13) (4.14) and (4.15) into (4.11) gives 

^(Vp,u) + ^||Vp|| 2 < CH| + Ce||Vti||. (4.16) 
at 2 

Multiplying (4.10) by D x suitably large and adding it to (4.16), we finally deduce the 
lemma since e > is sufficiently small. This completes the proof of Lemma 4.1. 

Our next goal is to deal the higher order estimate of (p, u). 

Lemma 4.2. For any < t < T, there exists a constant D 2 > sufficiently large 
which is independent of e, t and T such that 

f t {D 2 HMt),u(t)) + E <9?v P ,9?«)(t)} 

+0(1^^^)11? + l|v4t)Tl'lT< Ce(||Vp(t)|| 2 + ||^(t)|| 2 ), 
where Hi(p,u) is equivalent to ||V(p,w)(£)||!, if e is small enough. 

Proof: For each multi-index a with 1 < \a\ < 3, by applying d£ to (2.1)i — (2.1)2, 
multiplying them by <9"p, d"u respectively and then integrating them over M 3 , we have 

!!H^M|| 2 + a||9>|| 
= (%p, d*F) + <d>, d%G) (4.18) 

= h + h, 

where Jj, i — 1,2 are the corresponding terms in the above equation which will be esti- 
mated as follows. Here and in the sequel proof, the repeated index denotes summation 
over the index. 

Firstly, for Ji, it holds that, 



h = - (R+ * v v )Kl (d> d% (pV • u) ) - «i (<9> d£ (u • Vp) ) 

= - (fei^i (a^p, pd£ ( v • u ) ) - ki (a°p, ( u • v^p) ) 
- E" ° ^±^c£ (^p, ar Vf (v • «)> 



(4.19) 



Cv 

l/3|<l«|-l 

4 

= E J i,i- 
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By symmetry, Lemma 3.3 and some tedious but straightforward calculations, we can 
obtain 

X> < C||(p J « ja )|| 3 (||Vp||| + \\u\\l) < Ce(\\Vp\\l + \\u\\l). (4.20) 

i=2 

From (2.1)!, we have 

,. C v . _ (IK^OU 2 . 

div u = — — r {pt + mu ■ Vp H ), 

K 2 C V + (R + C v )Kip C v 

k 2 c v + (R + c v )Kip anjpu 2 

p t = (IIV U — K\U ■ \p , 

c v c v 

which imply that 

= (R + c v )(dZp,pdZ( £ (p t + Kl u ■ Wp + 



< 



^mpn, K2Cv+i £ +Cv)Kip ) + Ce(\\Vp\\l + IHIi) (4.21) 
< {J ^imP)\ K2Cv+i l Cv)KlP ) + Ce(\\Vp\\l + IHII). 

Combining (4.19) and (4.21), we have 

I < ^^<(^p) 2 , ; ) + C7e(||Vp||| + \\u\\l). 

2 at k 2 c v + {K + c v )Kip 

Taking the same argument to deal I 2 , then there exists a function Hi(p,u) which is 
equivalent to ||(Vp, Vn) j || and satisfies 

jHMt)Mt)) + llVtiHl < Cc(||Vp|||) + N| 2 ). (4.22) 

By a direct computation, we have the estimate on ||V<9"p|| for 1 < |a| < 2 as following 

y E ll Va ^H 2 + l E (d>^d: P )<Ce(\\Vp\\ 2 2 ) + \\u\\l). (4.23) 

l<a<2 l<|a|<2 

Since e is sufficiently small, multiplying (4.22) by D 2 suitably large and adding it to 
(4.23) give (4.17). Thus we completes the proof of the lemma. 

Finally, by symmetry and some tedious but straightforward calculation, we have 
the energy estimates on the entropy. 

Lemma 4.3. It holds that 

j t \\sm<C\\u(t)MW)\\l (4-24) 

for any < t < T. 
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5 The proof of global well-posedness 

In this section, we are devoted to prove Proposition 2.2. We first consider a priori 
decay-in-time estimates on (Vp, u). This will be based on Proposition 3.1 about the 
decay estimates on the linearized system. The decay-in-time estimate on p can be 
derived by decay-in-time estimates on (Vp, u). Precisely, we have the following lemma. 

Lemma 5.1. Let (p,u,s) be the solution of (2.1), then p satisfies the following in- 
equality 

||Vp|| < CK (1 + i)"t + Ct f\l + t - r)-I(||Vp(r)|| 2 + \\u(r)\\ 3 )dr, (5.1) 

Jo 

for any < t < T, where K = \ \(po,u )\\ L i nH 3 as in (2.5). 

Proof: From the Duhamel's principle and Proposition 3.1, we have 

||Vp|| <CK Q (l + t)-i+Ce f\l + t-T)-$(\\(F,G)(T)\\ L i nH i)dT. (5.2) 

Jo 

The nonlinear term (F, G) can be estimated as following 



(F,G)|| L i <C||(p,« )a )(t)||i(||Vp(0||i + |Kt)||2) 
<Ce(||Vp(t)|| 1 + |Kf)|| 2 ), 



(5.3) 
(5.4) 



IK^GOHi <C||(p ) « )a )(t)|| W ri.»(||Vp(0||2 + |Kt)||3) 
<C7e(||Vp(t)|| 2 + ||«(t)||3), 

Putting (5.3) and (5.4) into (5.2), we obtain (5.1). The proof of Lemma 5.1 is com- 
pleted. 

Now we are in a position to prove Proposition 2.2. 

Proof of Proposition 2.2. We do it by three steps. 

Step 1: Since e > is sufficiently small, from Lemma 4.1 and Lemma 4.2, we have 
a function H 2 (p,u) which is equivalent to ||(p, it) 1 1§ and satisfies 

j t H 2 (p(t)Mt)) + C(||Vp(f)||5 + |M| 2 3 ) < 0, (5.5) 
for any < t < T, which implies (2.3). 

Step 2: Multiplying (2.1) 2 by u, integrating over M 3 and using Cauchy-Schwarz 
inequality, we have 

^\\u\\ 2 + C\\u\\ 2 < C||Vp|| 2 + Ce||VM|| 2 . (5.6) 



dt 
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Now we define the temporal energy functional 

H 3 (t) = \\u(t)\\ 2 + D 2 H 1 (p(t),u(t))+ J2 <q?Vp,^u)(f) 

l<|a|<2 

for any < t < T, where it is noticed that H 3 (t) is equivalent to ||Vp||| + 
D2 can be large enough. 

Combining Lemma 4.2 and (5.6), we obtain 

jH 3 (t) + C{\\V 2 p(t)\\\ + IKOIIs) < C\\Vp\\ 2 . 
Adding ||Vp|| 2 to both sides of the above inequality gives 

j t H 3 (t) + D 3 H 3 (t)<C\\Vp\\ 2 , 
where D 3 is a positive constant independent of e. Set 

M(t) = sup (1 + t)§# 3 (t), 

0<T<t 

and notice that 

l|Vp(r)|| 2 + ||u|| 3 < Cy/H 3 (r) < C{\ + r)-*y/M(t), < r < t < T. 
Then it follows from Lemma 5.1 that 

||Vp(t)|| < CK (1 + t)-l + Ce | *(1 + t - r)-f (1 + r)-f dr^M{t) 
< C(l + t)- L l(K + e^M{t)). 

Hence, by the Gronwall's inequality, (5.7) and (5.10) lead to 

H 3 (t) < e- D ^H 3 (0) +C^e- D ^ t -^\\Wp(r)\\ 2 dr 

< e~ D ^H 3 (0) +Cf*e- D ^- T \l + T)'ldr(K 2 + e 2 M(t)) 

< C(l + t)-l(K 2 + e 2 M(t)). 

Since M(t) is non-decreasing, we have from (5.8) and (5.11) that 

M(t) < C(K 2 + e 2 M(t)), 
for any < t < T, which implies that 

M(t) < CK 2 , 

since e > is small enough. Thus we obtain (2.7) from (5.9) and (5.13). 
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Next, using Proposition 3.1, (5.3) and (5.4), it follows from the Duhamel's principle 
that 

||p(f)|| <CK (l + t)-* + Cj*(l + t-T)-^(F,G)(T)\\ LlnL 2dT 

< CK (1 + t)~l + C/ *(l + t - r)"i (||Vp(r)|| 2 + \\u\\ 3 )dr 

< CK (1 + t)~i + CK f*(l + t - r)-f (1 + t)~Ut 

<K (i + t)-i, 

for any < t < T. Thus (2.6) is proved. 

Step 3: By Lemma 4.3 and Gronwall's inequality, we arrive at 

s(t) < S (0)exp{C J*\\u(T)\\ 3 dr} 

<s(0)exp{CK J*(l + T)-Ur} 
< s(0) exp{CK } 

i.e., for any < t < T, 

\\(p,u,s)(t)\\l<C\\( Po ,u ,s )\\l exp{CK } . 

Hence (2.4) holds. For (2.8), using the above estimates and (2.1) we have 

\\d t (p,u,s)(t)\\ <C(|| M || 1 + ||Vp||) 
< CK (l + t)-l 

for any < t < T. Thus, (2.8) is proved and this completes the proof of Proposition 
2.2. 
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